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Summary 

Experimental  data  are  presented  for  the  development  of  compressible  turbulent  boundary  layers  moving 
from  a solid  onto  a transpired  surface,  and  vice  versa.  The  data  were  obtained  at  free  stream  Mach  numbers  of 
1-8  and  3-6  (Reynolds  numbers  per  metre  of  2-45  x 10^  and  545  X lO’  respectively)  and  transpiration  rates 
(pwVw/PiUi)  up  to  0-0044,  for  nominally  zero  heat-transfer  conditions. 

The  results  are  compared  with  theoretical  predictions  obtained  from  a method  based  on  the  eddy-viscosity 
concept. 

The  work  described  in  this  report  is  based  on  experimental  investigations  conducted  by  the  author  while  at 
Cambridge  University  Engineering  Department,  and  presented  in  a Ph.D.  dissertation  (Marriott*). 
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1.  Introdnctioa 

Papers  by  Jeromin^  and  Squire^  have  reported  experiments  concerning  the  development  of  compressible 
turbulent  lx>undary  layers  on  surfaces  with  uniform  transpiration.  A logical  extension  of  this  work  is  the 
investigation  of  the  situation  where  transpiration  rates  are  non-uniform.  The  simplest  case  of  non-uniform 
transpiration,  namely  a plane  discontinuity,  where  a solid  surface  boundary  layer  encounters,  and  develops  on, 
a transpired  surface,  is  the  subject  of  the  present  investigation. 

Apart  from  the  work  of  Kays  and  others  at  Stanford  (summarised  by  Kays*),  and  that  of  McQuaid^,  studies 
of  non-uniform  transpiration  are  generally  insubstantial. 

The  purpose  of  the  pre^nt  investigation  was  to  obtain  basic  experimental  data  for  a compressible  turbulent 
boundary  layer  developing  over  a straight  discontinuity  in  surface  transpiration  rate,  for  which  tabulated 
profile  data  are  presented  in  this  report.  A theoretical  eddy  viscosity  method,  due  to  Verma^,  was  used  for 
predictions  of  the  flow  developments,  and  these  are  also  presented  in  comparison  with  the  experimental 
results. 

The  eddy-viscosity  formulation  is  given  in  Appendix  A and  the  parameters  used  in  the  method  are  given  in 
Table  1. 


2.  Experimental  Details 

2.1.  Apparatus 

All  the  experimental  work  in  the  present  study  has  been  carried  out  in  an  intermittent  blow-down 
supersonic  tunnel  at  the  Cambridge  University  Engineering  Department.  A description  of  the  tunnel  and 
associated  transpiration  apparatus  has  been  given  by  Jeromin^  and  Squire^.  A photograph  of  the  tunnel  can  be 
seen  in  Fig.  1,  and  a sketch  of  the  working  section  in  Fig.  2. 

At  the  two  Mach  numbers  investigated,  Mach  1*8  and  3-6,  the  stagnation  chamber  pressures  were 
respectively  82*7  kN/m^  and  887  kN/m^  gauge.  There  was  no  provision  to  control  the  stagnation  tempera- 
ture, which  was  found  to  vary  by  about  ±2  K during  a run,  having  a nominal  value  of  293  K±4K.  The 
temperature  of  the  transpired  air  was  close  to  stagnation  temperature  or  slightly  below. 

The  working  section  of  the  tunnel,  as  described  in  Jeromin^,  allows  for  boundary-layer  measurements  to  be 
made  over  a surface  approximately  200  nun  long  by  100  mm  wide.  This  surface  forms  part  of  a plate, 
measuring  approximately  400  nun  x 100  mm,  which  can  be  easily  interchanged,  allowing  for  a variety  of 
surface  configurations  to  be  tested. 

The  porous  surface  in  the  present  investigation  was  constructed  from  rolled  sintered  bronze,  butt-jointed  to 
the  solid  portion  of  the  plate,  as  shown  in  Figs.  3 and  4.  Copper-constantan  type  thermocouples  were  inserted 
in  the  plate  for  the  measurement  of  wall  temperature.  For  the  porous  to  solid  measurements,  the  plate  was 
simply  reversed.  The  zero  transpiration  data  were  obtained  on  a completely  solid  surface. 

2.2.  Data  Measurement 

Boundary-layer  total-pressure  measurements  were  obtain  using  a flattened  pitot  probe,  as  shown  in  Figs.  5 
and  6,  whidi  was  attached  to  a pressure  transducer  and  mounted  on  a traversing  mechanism.  A sliding 
potentiometer  on  the  traverse  mechanism  provided  an  electrical  signal  which  was  fed  into  one  channel  of  an 
X~  Y recorder,  the  other  channel  being  connected  to  the  pressure  transducer,  thus  enabling  a complete  profile 
record  to  be  obtained  during  a single  tunnel  run. 

Temperature  profiles  were  obtained  in  a similar  manner  using  a copper-constantan  thermocouple  embed- 
ded in  a 10  degree  included-angle  steel  cone,  itself  mounted  on  ceramic  tube  to  reduce  heat-transfer  effects, 
and  connected  directly  to  the  X~Y  recorder.  Dimensions  can  be  found  in  Fig.  6. 

13.  Data  Reduction 

All  the  X-Y  recorder  plots  were  analysed  by  reading  off  a series  of  points,  each  corresponding  to  a 
measured  pressure  (or  temperature)  and  its  associated  distance  from  the  surface,  together  with  the  relevant 
scale  lengths.  These  results  were  fed  into  a computer  program  which  calculated  the  velocity  and  static 
temperature  profiles  from  the  normal  Rayleigh  pitot  formulae,  using  static  pressure  from  static  probe 
measurements.  About  SO  points  were  found  adequate  to  cover  the  complete  profile,  and  these  were  spaced 
according  to  the  rate  of  change  of  pressure  (or  temperature),  more  being  taken  in  the  region  near  the  wall, 
where  the  rate  of  change  was  normally  largest. 


3 


From  these  basic  results,  boundary-layer  parameters,  such  as  displacement  thickness  and  momentum 
thickness,  were  calculated.  Shear-stress  profiles  were  also  obtained,  using  the  momentum  equation  and  the 
measured  velocity  and  temperature  profiles. 

2.4.  Two-Dimensioiiality 

Jeromin^,  who  used  essentially  the  same  equipment  as  used  in  the  present  tests,  measured  velocity  profiles 
25  mm  either  side  of  the  centre  line,  and  reported  variations  of  6 and  5*  of  about  2 per  cent,  from  which  he 
assumed  that  three-dimensional  effects  were  negligible.  However,  in  a further  publication,  Jeromin^  con- 
cluded there  may  well  be  three-dimensional  effects  present  at  Mach  3*6,  since  he  found  that  skin-friction 
coefiScients  evaluated  from  the  momentum-integral  equation  were  about  8 per  cent  too  high  when  compared 
to  the  results  of  other  investigators.  His  skin-friction  results  evaluated  at  Mach  1*8  and  2*5  agreed  well  with 
published  values,  suggesting  that  in  these  cases,  three-dimensional  effects,  if  any,  were  minimal. 

In  the  present  investigation,  spanwise  pressures  were  measured  simultaneously  using  a number  of  pitot 
probes  and  typical  results  are  shown  in  Fig.  7.  Compared  with  the  results  at  Mach  1*8,  those  at  Mach  3*6  show 
significant  total  pressure  variations  across  the  span.  As  will  be  shown  later,  the  experimental  Mach  3*6  data  are 
in  good  agreement  with  theory,  and  it  is  therefore  considered  that  the  error  introduced  by  three-dimensional 
effects  is  not  large. 

2.5.  Test  Condidons 


Solid-Porous 

Porous-Solid 

Nominal  Mach  number 
Stagnation  pressure 

Stagnation  temperature 
Reynolds  number  per  metre 

1*8 

82*7  kN/m" 
gauge 
293  K 
2*45x10^ 

3*6 

887kSlm^ 
gauge 
293  K 
5*45x10^ 

3*6 

887  kN/m^ 
gauge 
293  K 
5*45  xlO"" 

Transpiration  rate  — — 
P\V\ 

0 

0 

0*00174 

0*00311 

0*00446 

0*00164 

0*00294 

0*00421 

0*00294 

With  the  discontinuity  as  the  datum  the  traverse  positions  were: 


-25*4  mm. 

0 mm, 

3*2  mm. 

6*3  mm. 

9*5  mm. 

12*7  mm. 

19*0  mm. 

25*4  mm. 

38*1  mm. 

50*8  mm. 

63*5  mm. 

88*9  mm. 

114*3  mm. 

139*7  mm. 

171*4  mm. 

3.  Discussion  of  Results 
3.1.  Free-Stream  Mach  Number  Distribution 

Although  the  experimental  results  were  obtained  in  a nominally  zero  pressure  gradient,  some  localised 
pressure  gradients  were  present,  as  is  apparent  from  the  free-stream  Mach-number  distributions  shown  in 
Figs.  8 and  9.  These  were  caused  by  the  shock  fields  in  the  working  section,  which  can  be  seen  in  the  schlieren 
photographs  in  Fig.  5 and  also  by  the  presence  of  three-dimensional  effects  described  previously.  The  most 
serious  flow  disturbance  was  generated  by  the  discontinuity  itself,  in  that  the  sudden  change  in  boundary-layer 
displacement  thickness  downstream  of  the  discontinuity  gives  rise  to  an  expansion  wave  (porous  to  solid 
change),  or  a compression  wave  (solid  to  porous).  This  can  be  clearly  seen  in  the  schlieren  photographs,  and  is 
apparent  in  the  results  by  a sudden  change  in  free-stream  Mach  number.  The  fall  in  free-stream  Mach  number 
for  the  Mach  1 *8  flows,  at  about  125  mm  downstream  of  the  discontinuity,  is  due  to  a reflected  shock  starting 
from  the  leading  edge  of  the  plate,  and  visible  in  the  schlieren  pictures. 
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Due  to  the  decrease  in  shock  angle  with  increasing  Mach  number,  the  effects  of  the  discontinuity-induced 
shock  wave  are  felt  further  downstream  for  the  flows  at  Mach  3*6  than  those  at  Mach  1*8.  On  average,  the 
affected  regions  extend  to  about  12  mm  and  25  nun  for  the  Mach  1'8  and  3-6  flows  respectively. 


3.2.  Veiodty-Profile  Developments 

Figs.  10  to  20  show  typical  profile  developments,  together  with  theoretical  predictions  obtained  using  an 
eddy- viscosity  method  due  to  Verma*.  The  experimental  velocity  and  temperature  profiles  at  x = -25*4  mm 
are  used  in  the  method  to  provide  starting  conditions  for  the  theoretical  developments.  Experimentally 
measured  free-stream  temperature  and  Mach-number  distributions  are  used  in  the  method  for  the  free-stream 
boundary  conditions.  It  can  be  seen  that  the  experimental  and  theoretical  results  agree  very  well,  particularly 
at  the  lower  transpiration  rates.  The  major  differences  to  be  found  between  the  experimental  and  theoretical 
profiles  occur  at  the  rear  of  the  plate,  most  noticeably  at  the  higher  transpiration  rates.  This  is  apparent  in  Fig. 
18,  where  for  the  Mach  3-6,  F= 0-00421  development,  the  measured  velocity  profile  exhibits  an  inflexion 
typical  of  that  in  a boundary  layer  near  to  separation.  Indeed,  it  was  found  that  separation,  or  ‘blow-off’, 
conditions  did  exist  at  the  last  measuring  station  (at  x = 171-4  mm),  the  results  for  which  are  not  shown,  for 
both  the  F=  0-00421  and  F= 0-00294  flows  at  Mach  3-6. 

Some  individual  experimental  profiles  are  also  seen  to  be  in  slight  disagreement  with  the  theoretical 
predictions.  This  is  possibly  a result  of  localised  shock-wave  interference  at  the  outer  edge  of  the  boundary 
layer,  introducing  errors  in  the  determination  of  the  free-stream  conditions,  and  hence  producing  unrealistic 
profile  shapes. 

It  will  also  be  observed  that  the  generally  good  agreement  between  experiment  and  theory  for  the  Mach  3-6 
developments  suggests  that  the  three-dimensional  effects  mentioned  earlier  are  minimal. 

3.3.  Temperatwe-Profile  Developments 

Figs.  11  to  21  also  show  typical  temperature  profile  develoments,  together  with  theoretical  predictions 
(using  the  method  due  to  Verma^).  Good  agreement  is  obtained  between  the  predicted  and  experimental 
profiles,  but  deteriorating  towards  the  rear  of  the  plate,  particularly  at  the  higher  blowing  rates.  Again,  as  was 
discussed  in  the  previous  section,  impending  ‘blow-off’  is  probably  responsible  for  the  inflexion  in  the  profile 
shapes.  It  is  interesting  to  observe,  however,  that  the  inflexion  is  also  apparent  in  the  predicted  profiles. 

Figs.  22  and  23  show  typical  distributions  ofT/Ti  plotted  against  u/Ui,  together  with  the  Crocco  relation,  in 
which  the  recovery  temperature,  T„  was  assumed  equal  to  the  wall  temperature  Th,.  Thus  it  can  be  written; 


T uy 

r,  Ti  \ Ti  )\uj 


(1) 


While  the  agreement  between  the  results  and  the  Crocco  relation  is  quite  good  at  Mach  1 -8,  at  Mach  3-6  it  is 
less  satisfactory.  It  should  be  noticed,  however,  that  these  comparisons  are  somewhat  biased,  in  that  the 
temperature  profiles  in  the  region  0<u/(/i<0-7  approximately,  have  been  deduced  from  a simple  linear 
extrapolation  of  the  recorded  probe  temperature  nearest  the  wall,  to  that  of  the  wall  itself. 


3.4.  Integral  Parameters 

In  general,  the  predicted  developments  of  0 compare  well  with  the  experimental  results,  as  can  be  seen  in 
Figs.  24  to  26.  It  might  be  expected  that  some  disagreement  would  exist  just  downstream  of  the  discontinuity 
where  there  are  marked  changes  in  the  values  of  0,  and  which  can  be  attributed  to  the  interference  effects  on 
the  boundary-layer  of  the  shock  wave  produced  at  the  start  of  the  transpired  region.  The  shock-wave 
boundary-layer  interaction  may  also  account  for  the  slight  differences  noticeable  between  the  experimental 
and  predicted  developments,  since  modifications  to  the  initial  structure  of  the  boundary-layer  will  obviously 
affect  the  remainder  of  its  growth. 

The  results  for  the  porous-to-solid  situation  show  most  disagreement  with  theory,  as  do  the  velocity  and 
temperature  profile  predictions  for  this  particular  case.  It  can  be  seen  that  the  positions  of  shock-wave 
interferences  on  the  boundary  layer,  shown  in  the  Schlieren  photographs,  correspond  closely  to  those  where 
there  are  marked  changes  in  the  developments  of  0.  It  therefore  appears  probable  that,  as  before,  shock-wave 
boundary-layer  interactions  are  responsible. 

Figs.  27  and  28  show  composite  plots  of  shape  factor  variations  against  for  the  various  Mach  number 
and  transpiration  rates,  together  with  some  equilibrium  boundary-layer  results  of  Squire*.  (For  clarity. 
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individual  experimental  points  have  been  omitted  and  replaced  by  mean  curves.)  The  results  are  consistent 
with  those  of  Squire,  except  for  the  solid  surface  Mach  3-6  development,  which  show  slightly  different  trends. 
This  might  possibly  be  due  to  three-dimensional  effects,  or  to  the  slightly  different  flow  conditions  for  the  two 
sets  of  experiments.  (In  the  present  investigation,  the  duralumin  liner  forming  the  nozzle  in  the  tunnel  was 
moved  152  mm  further  upstream,  relative  to  its  position  in  Squire’s  experiments,  allowing  a longer  flow 
development.  It  is  possible,  therefore,  that  there  existed  small  differences  in  the  flow  conditions  for  the  two 
configurations.) 

The  data  suggest  that,  in  the  present  experiment,  equilibrium  conditions  have  not  been  reached,  in  the  sense 
that  H has  not  become  a unique  function  of  Rt  and  F,  even  at  the  last  measuring  station,  since  the  curves  do  not 
attain  the  levels  inferred  from  Squire’s  results.  This  conclusion  is  supported  by  the  development  of  the 
y/$  vs.  u/Ui  profiles,  shown  in  I^gs.  29  to  31,  in  which  it  can  be  seen  that  the  experimental  profiles  do  not  reach 
a common  shape  at  the  end  of  the  development.  However,  it  would  appear  that  the  flow  is  in  a near- 
equilibrium condition,  since  the  profiles  show  a reasonable  degree  of  collapse  at  the  rear  of  the  plate  for  the 
solid-to-porous  flows,  and  the  one  comparison  available  of  an  equilibrium  profile  from  Squire’s  results 
confirms  this. 

Anomalously,  the  y/$  vs.  u/Ui  curves  for  the  porous-to-solid  configuration  tend  to  suggest  that  an 
equilibrium  layer  is  established  in  a relatively  short  distance  downstream  of  the  discontinuity,  even  though  the 
‘equilibrium’  profile  shape  does  not  agree  with  either  the  solid-surface  equilibrium  profile  of  Squire,  or  the 
soUd-surface  results  from  the  present  investigation.  This  may  be  due  to  enors  in  the  experimental  d variations 
over  the  rear  of  the  plate,  which  appear  to  be  too  high  compared  to  the  theory.  Certainly,  for  this  case,  the  Hvs. 
R»  curves  suggest  that  equilibrium  conditions  have  yet  to  be  reached. 

3^.  Skfai  nMioa  Variattons 

I^gs.  32  and  34  show  the  predicted  skin  friction  developments,  while  Figs.  33  and  35  show  the  variation  in 
dynamic  pressure  measured  by  positioning  the  total-pressure  probe  onto  the  surface  as  a Stanton  tube.  It  is 
interesting  to  note  that  the  variations  in  both  predicted  skin  friction,  and  measured  dynamic  pressure  near  the 
surface,  correspond  remarkably  well,  and  indicate  that  most  of  their  variations  occur  within  the  first  25  mm 
downstream  of  the  discontinuity  (about  4 to  5 boundary-layer  thicknesses).  The  exception  is  the  porous-to- 
solid  development,  which  appears  to  vary  rather  more  slowly,  but  again,  the  trend  predicted  by  the  theory,  and 
that  shown  by  the  dynamic-pressure  variation,  are  in  good  agreement. 

3.6.  Inner-Law  Profile  Representation 

Using  the  inner  law  proposed  by  Squire^,  of  the  form 


+B 


(2) 


where 


G{u 


V[F«*+(Q/2)] 


du\ 


U,' 


'■-7. 


and  the  skin-friction  variation  given  by  theory,  inner-law  plots  were  obtained,  as  shown  in  Figs.  36  to  41.  It  can 
be  seen  that  the  profiles  (with  the  exception  of  the  porous-to-solid  and  the  M = 3*6,  F = 0-00421  cases)  show 
an  extremely  good  collapse,  bearing  in  mind  the  uncertainty  in  value  of  Q.  Dunbar*  reported  that  a variation 
of  10  per  cent  in  the  value  of  Q can  give  rise  to  an  overall  change  in  the  levels  of  G(u' ) of  about  5 per  cent,  and 
thus  errors  in  Q could  well  account  for  the  small  variations  between  profiles,  particularly  just  downstream  of 
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the  discontinuity  where  rapid  changes  occur.  The  poor  result  at  Af  = 3*6,  F = 0*0042 1 could  also  be  due  to 
large  errors  in  Q,  which  for  this  case  is  small  (approaching  zero  at  the  separation  point). 

Similar  arguments  could  also  explain  the  lack  of  collapse  for  the  porous-to-solid  profiles.  In  this  case  it  will 
be  noticed  that  the  maximum  levels  of  G(u^)  for  the  profiles  in  the  region  of  the  discontinuity  are  far  higher 
than  those  obtained  at  the  rear  of  the  solid-to-porous  configiuation  with  the  same  transpiration  rate. 

With  the  exception  of  the  porous-to-solid  configuration,  the  results  for  the  solid-to-porous  cases  show  that 
the  inner  part  of  the  boundary  layer  responds  very  quickly  to  the  change  in  surface  conditions. 

3.7.  Shear-Stress  Profile  Calculations 

A numerical  technique  was  used  for  calculating  the  shear-stress  distributions  from  the  measured  velocity 
and  temperature  profiles,  using  the  momentum  equation.  The  momentum  equation,  usually  written 

T(y)-T«,-y^  = j^  ^(pu^)dy'-u(y)  -^(pu)  dy' + p^V^u(y),  (3) 

is  rewritten  in  the  form 

T(y)-T»,-y^=^^  (pM*)</y'-M(y)  (pw)<fy'+|^  ||^  (pu) dy'j  </«'j+p»,VXy).  (4) 


Whilst  the  two  equations  are  essentially  the  same,  the  computational  advantage  in  the  second  lies  in  the 
single  numerical  difierentiation. 

Using  equation  (4),  the  shear-stress  profiles  were  calculated,  and  typical  cases  are  shown  in  Figs.  42  to  44.  It 
can  be  seen  that  the  experimental  results  agree  reasonably  well  with  the  overall  predictions,  considering  the 
limitations  of  the  calculation  technique.  Shock-wave  boundary-layer  interactions  producing  profile  distortions 
obviously  introduce  errors,  and  is  particularly  relevant  just  downstream  of  the  discontinuity.  In  addition, 
errors  in  the  skin-friction  value  used  are  reflected  directly  in  the  shear-stress  profiles. 

Analytical  errors  are  introduced  primarily  in  obtaining  the  x -derivatives  in  the  calculation  procedure.  To 
reduce  these  errors,  the  derivatives  were  found  by  fitting  a curve  to  the  function  inside  the  square  brackets  in 
equation  (4).  However,  the  errors  may  become  dominant  towards  the  edge  of  the  layer,  where  small  errors  in 
the  shape  of  the  profiles  give  rise  to  finite  shear-stress  values.  A non-zero  shear  stress  at  the  outer  edge  of  the 
boundary  layer  was  assumed,  for  the  purposes  of  the  numerical  procedure,  to  be  a result  of  a localised  pressure 
gradient,  which  was  then  applied  proportionately  throughout  the  layer.  In  this  manner,  the  shear-stress 
profiles  were  forced  to  zero  at  y = S{u/U\  - 1). 

Predicted  and  experimental  shear-stress  profiles  for  the  porous-solid  situation,  although  exhibiting  similar 
shape,  are  not  in  good  agreement  at  the  start  of  the  development.  However,  the  agreement  improves  towards 
the  end  of  the  plate  where  it  is,  in  fact,  extremely  good.  This  behaviour  might  be  expected  following  the  rather 
poor  results  shown  in  the  logarithmic  inner-law  plots,  with  the  associated  uncertainty  over  the  magnitude  of 
the  skin  friction. 


3A  IVfixkag-Lcagth  DistribotioiU 

Using  the  shear-stress  profiles  calculated  as  described  in  the  previous  section,  the  mixing  length  distribU' 
tions  were  obtained,  from 


'-O'/t 


(5) 


and  shown  in  Figs.  45  to  49.  In  general,  the  results  show  a reasonable  collapse  onto  a straight  line  in  the  wall 
region,  with  a slope  of  about  0*4.  However,  at  the  higher  transpiration  rates,  and  particularly  for  the  mixing 
length  dutributions  calculated  from  profiles  in  the  region  just  downsteam  of  the  discontinuity,  some  scatter  in 
the  results  is  evident.  This  is  to  be  expected,  since  the  shear-stress  profiles  in  this  region  showed  some 
disagreenwnt  with  the  theoretical  predictions,  particularly  near  the  wall. 

The  porous-to-solid  results  do  not  show  a good  collapse,  which  is  to  be  expected  considering  the 
dissimilarity  between  the  theoretical  and  calculated  shear-stress  profiles. 
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In  the  outer  region,  the  mixing-length  distributions  do  not  show  any  consistent  trends.  The  scatter  can  be 
attributed  to  errors  in  both  the  shear-stress  profile  and  the  slope  of  the  velocity  profile  (du/dy)  in  the  outer 
region  (particularly  just  downstream  of  the  discontinuity  where  shock-wave  interaction  effects  are  dominant), 
both  parameters  tending  to  zero  at  the  outer  edge  of  the  layer.  However,  results  suggest  that  the  adoption  of  a 
constant  value  for  //£  if  between  0-07  and  0-08  is  not  uru'easonable. 

1 he  outer  region  mixing-length  distributions  for  the  porous-to-solid  case  appear  to  be  somewhat  lower  than 
those  of  the  solid-to-porous  developments.  However  this  might  be  expected,  considering  the  experimental 
shear-stress  profiles  for  this  case  which  also  appear  low. 

3.9.  Eddy-Viscosity  Distributioiis 

These  were  obtained  from  the  relation 


T du 
»»,  = - — 
P dy 


(6) 


and  shown  in  Figs.  45  to  49.  The  profiles  exhibit  a fair  collapse  in  the  outer  region,  except  for  the  first  few 
profiles  immediately  downstream  of  the  discontinuity.  Apart  from  these,  a value  of  about  O-OIS  would  appear 
to  be  a reasonable  estimate  of  the  maximum  value  of  the  eddy  viscosity.  The  exceptions  are  the  results  for  the 
porous-to-solid  development,  which  appear  to  be  extremely  low,  having  a maximum  of  about  0-08.  This 
particular  case  reflects  the  experimental  shear-stress  profiles,  which  can  be  seen  to  be  much  lower  than  those 
given  by  prediction.  From  these  findings,  it  would  seem  probable  that  there  are  deficiencies  in  the  outer  region 
eddy-viscosity  model  as  used  to  predict  the  porous-to-solid  development. 


4.  Condusioiis 

4.1.  The  boundary  layer  was  not  found  to  reach  an  equilibrium  state  even  at  the  end  of  the  measuring  region 
(equivalent  to  about  35  boundary-layer  thickness).  However,  from  the  experimental  results  for  the  Mach  1-8 
flows  at  the  lower  transpiration  rates,  it  would  appear  that  a near  equilibrium  condition  is  attained,  as  indicated 
by  the  near-collapse  of  the  u/Ui  vs.  yjd  curves. 

4.2.  Except  for  the  porous-to-solid  situation,  the  inner  part  of  the  boundary  layer  appears  to  respond  rapidly 
to  the  change  in  surface  conditions,  as  is  shown  by  the  generally  good  collapse  of  the  logarithmic  law  in  this 
region.  Lack  of  collapse  for  the  porous-to-solid  development  is  not  necessarily  indicative  of  a slow  response  of 
the  inner  region.  Errors  in  the  magnitude  of  the  skin  friction  may  be  responsible  for  this. 

4.3.  Uncertainties  in  the  experimental  results,  particularly  for  the  value  of  skin  friction,  makes  quantitative 
analysis  of  the  mixing-length  and  eddy-viscosity  results  difficult.  However,  it  would  appear  that  the  adoption 
of  a mixing-length  relation  of  the  form  I = const,  x y in  the  inner  region  is  not  unreasonable.  In  addition,  while 
the  eddy-viscosity  results  do  not  show  a good  collapse  in  the  outer  part  of  the  flow,  there  would  appear  to  be  no 
justification  in  adopting  a complex  model  in  this  region. 

4.4.  In  general,  the  experimental  results  agree  reasonably  well  with  the  eddy  viscosity  theoretical  predictions. 
Most  noticeable  differences  occur  in  the  regions  of  shock  wave  boundary  interactions,  and  towards  the  end 
of  the  measuring  region  where  the  boundary  layer  was  tending  towards  separation,  particularly  at  the  higher 
transpiration  rates. 

4.5.  Three-dimensional  effects  may  be  present  in  the  Mach  3*6  flows,  although  their  influence  on  the  data  is 
considered  minimal. 

5.  Tables  of  Experimental  Results 

5.1.  WaO  Temperatnre 

The  three  thermocouples  in  the  porous  surface  each  gave  approximately  the  same  reading,  being  nominally 
2 K below  stagnation  temperature  for  the  Mach  1-8  flows,  and  6 K below  at  Mach  3*6. 
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5.2.  fV«e-Slream  Velodty  and  Temperature 

Tables  2 and  3 show  the  measured  free-stream  velocities  and  temperatures.  Identification  of  the  flow 
situation  is  given  at  the  head  of  each  column,  in  the  form 

A-B-C 

where  A is  the  nominal  free-stream  Mach  number 
B is  the  flow  situation,  set  equal  to 

0 for  the  solid  surface 

1 for  the  solid-to-porous  surface 

2 for  the  porous-to-solid  surface 

C is  the  transpiration  rate,  F (=PwVw/p\U\) 

5.3.  Displacement  and  Momentum  Thicknesses 

Tables  4 and  5 give  the  displacement  and  momentum  thicknesses  respectively,  calculated  from  the' velocity 
and  temperature  profiles. 

5.4.  Tabulated  Velocity  and  Temperature  Profiles 

The  remaining  tables  list  the  values  otu/Ui  and  T/  T i which  have  been  interpolated  onto  a fixed  y -grid  from 
the  experimental  results.  Each  set  of  results  is  again  identified  in  the  manner  described  in  the  previous  section. 

It  should  be  noted  that  the  data  points  shown  in  the  figures  are  those  from  experimental  and  not  the 
interpolated  values  given  in  the  tables. 
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UST  OF  SYMBOLS 


Ao 

B 

C, 

F 

H 

*l.*2 

/ 

M 

P 

Pr 

Pr. 

T 

T, 

u,U 

Ur 

V 

X 

y 

s 

8* 

St 

e 

p 


Damping-function  constant 

Log-law  constant 

Skin  friction  { = tJp\ U\) 

Transpiration  parameter  (=p„VJpiUi) 
Shape  factor 

Constants  in  eddy-viscosity  model 

Mixing  length 

Mach  number 

Static  pressure 

Laminar  Prandtl  number 

Turbulent  Prandtl  number 

Static  temperatiu-e 

Recovery  temperature 

Streamwise  velocity 

Friction  velocity  ( = y/rjpw) 

Velocity  normal  to  surface 
Streamwise  coordinate 


Normal  coordinate 

Boundary-layer  thickness  (at  u/Ui  = 0-995) 
Displacement  thicknesss  ^ = | (l  - ^y) 

Incompressible-displacement  thickness  ( = 1 


Momentum  thickness 


Kinematic  viscosity 


v. 

T 

P 


e 


Subscripts 

w 

1 

a 


‘Effective’  kinematic 


Shear  stress 
Density 


Value  at  wall 
Value  at  free  stream 
Atmospheric  conditions 


TABLE  1 


Parameters  Used  in  the  Eddy- Viscosity  Method 

Ao  = 26  (constant  in  Van  Driest  damping  function) 
ki  = 0-4  (inner-region  constant) 
jtj  = 0-016  (outer-region  constant) 

Pr  = 0‘7  (laminar  Prandtl  number) 

Pr,  = 0-89  (turbulent  Prandtl  number) 


TABLE  2 

Tables  of  Ut  (m/s)  at  Experimental  Stations 


X \ 
nun  \ 

1-8-0- 

0-0000 

1-8-1- 

0-00174 

1-8-1- 

0-00311 

1-8-1- 

0-00446 

3-6-0- 

0-0000 

3-6-1- 

0-00164 

3-6-1- 

0-00294 

3-6-1- 

0-00421 

3-6-2- 

0-00294 

-25-4 

— 

481-22 

659-05 

655-68 

659-82 

655-71 

657-70 

0 

482-56 

482-56 

484-96 

656-76 

656-59 

652-99 

658-52 

3-2 

655-74 

652-63 

653-85 

658-65 

6-3 

486-47 

484-48 

655.35 

654-34 

655-74 

659-66 

9-5 

482-31 

482-87 

476-69 

651-25 

653-68 

657-96 

651-74 

12-7 

481-25 

476-85 

655-65 

656-23 

656-92 

657-93 

653-56 

190 

482-31 

479-23 

477-56 

659-30 

654-47 

659-30 

656-62 

25-4 

475-86 

476-20 

474-52 

652-12 

654-63 

647-00 

657-96 

38- 1 

481-22 

482-59 

655-06 

652-65 

651-46 

660-91 

50-8 

483-15 

484-11 

653-77 

653-57 

652-85 

659-35 

63-5 

489-08 

484-94 

483-92 

483-92 

656-89 

657-53 

648-64 

651-38 

662-76 

88-9 

489-33 

478-80 

483-51 

655-62 

654-68 

651-76 

655-53 

114-3 

475-58 

487-42 

483-24 

482-70 

658-23 

654-18 

655-17 

651-14 

658-02 

480-15 

481-65 

479-91 

654-86 

653-97 

650-29 

661-80 

478-15 

474-57 

465-82 

467-47 

659-76 

656-46 

661-32 
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TABLE  3 


Tables  of  Ti  (Kelvin)  at  Experimental  Stations 


X \ 

mm  \ 

1-8-0- 

0-0000 

1-8-1- 

0-00174 

1-8-1- 

0-00311 

1-8-1- 

0-00446 

3-6-0- 

0-0000 

3-6-1- 

0-00164 

3-6-1- 

0-00294 

3-6-1- 

0-00421 

3-6-2- 

0-00294 

-25-4 

176-6 

WSSi 

83-0 

82-1 

83-7 

82-7 

84-3 

0 

mSm 

83-2 

83-3 

82-5 

83-8 

3-2 

177-1 

WSm 

83-2 

82-1 

82-6 

83-8 

6-3 

176-7 

mam 

82-8 

82-8 

83-2 

83-8 

9-5 

179-5 

WSlm 

81-8 

82-7 

83-6 

81-5 

12-7 

175-4 

175-2 

Ittll 

179-8 

82-7 

83-4 

83-6 

83-6 

83-3 

190 

179-5 

■Bl 

182-1 

84-0 

82-2 

84-0 

84-3 

25-4 

177-1 

178-6 

181-0 

84-6 

87-1 

86-9 

84-2 

38- 1 

177-3 

177-7 

178-8 

84-1 

85-4 

87-2 

81-9 

50-8 

174-7 

177-6 

177-4 

84-8 

86-6 

88-1 

81-2 

63-5 

175-6 

177-3 

178-8 

178-8 

83-7 

86-1 

85-4 

87-8 

82-1 

88-9 

178-2 

170-9 

175-2 

84-3 

85-4 

86-5 

80-6 

114-3 

180-6 

177-9 

176-3 

178-2 

83-4 

83-5 

85-7 

87-4 

81-9 

139-7 

179-8 

179-7 

182-1 

83-7 

85-3 

88-9 

83-1 

171-4 

180-1 

179-3 

183-9 

187-7 

83-7 

83-9 

83-9 

TABLE  4 


Tables  of  S*(nim)  at  Experimental  Stations 


1-8-0- 

1-8-1- 

1-8-1- 

1-8-1- 

3-6-1- 

3-6-1- 

3-6-1- 

1 

1 

3-6-2- 

0-0000 

0-00174 

0-00311 

0-00446 

0-0000 

0-00164 

0-00294 

0-00421 

0-00294 

-25-4 

1-209 

1-219 

1-255 

1-232 

2-736 

2-720 

2-705 

2-703 

7-361 

0 

1-288 

1-303 

1-295 

2-797 

2-797 

2-786 

7-607 

3-2 

1-293 

1-318 

1-321 

2-858 

2-858 

2-827 

7-587 

6-3 

1-290 

1-341 

1-351 

2-819 

2-835 

2-835 

7-638 

9-5 

1-377 

1-384 

1-488 

2-812 

2-807 

2-817 

7-597 

12-7 

1-308 

1-415 

1-488 

1-565 

2-883 

2-819 

2-822 

2-868 

7-899 

19-0 

1-478 

1-575 

1-643 

2-908 

3-025 

3-061 

7-940 

25-4 

1-534 

1-638 

1-775 

3-096 

3-299 

3-553 

7-922 

38-1 

1-641 

1-753 

1-877 

3-287 

3-612 

3-993 

7-297 

50-8 

1-694 

1-877 

2-106 

3-498 

3-957 

4-417 

7-224 

63-5 

1-400 

1-791 

2-045 

2-294 

3-099 

3-658 

4-282 

4-874 

7-219 

88-9 

1-971 

2-342 

2-715 

3-978 

4-793 

5-677 

7-448 

114-3 

1-737 

2-195 

2-611 

3-124 

3-261 

4-397 

5-380 

6-751 

7-470 

139-7 

2-410 

2-934 

3-635 

4-722 

6-043 

7-943 

7-777 

171-4 

1-750 

2-784 

3-813 

5-674 

3-490 

5-593 

7-803 
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TABLE  5 


Tables  of  0(min)  at  Experimental  Stations 


A 

mm  \ 

18-0- 

0-0000 

1-8-1- 

0-00174 

1-8-1- 

0-00311 

1-8-1- 

0-00446 

3-6-1- 

0-0000 

3-6-1- 

0-00164 

3-6-1- 

0-00294 

3-6-1- 

0-00421 

3-6-2- 

0-00294 

-25-4 

0-424 

0-427 

0-439 

0-427 

0-401 

0-381 

0-396 

0-389 

0-869 

0 

0-450 

0-455 

0-465 

0-404 

0-404 

0-401 

0-947 

3-2 

0-460 

0-452 

0-455 

0-411 

0-396 

0-399 

0-953 

6-3 

0-452 

0-470 

0-455 

0-394 

0-378 

0-389 

0-975 

9-5 

0-493 

0-480 

0-521 

0-396 

0-394 

0-384 

0-963 

12-7 

0-460 

0-485 

0-523 

0-536 

0-422 

0-399 

0-394 

0-386 

1-021 

190 

0-523 

0-539 

0-561 

0-411 

0-404 

0-419 

1-069 

25-4 

0-526 

0-554 

0-592 

0-445 

0-472 

0-493 

1-069 

38- 1 

0-564 

0-594 

0-627 

0-470 

0-505 

0-551 

0-975 

50-8 

0-574 

0-638 

0-681 

0-500 

0-551 

0-605 

0-950 

63-5 

0-500 

0-597 

0-671 

0-726 

0-462 

0-523 

0-582 

0-653 

0-958 

88-9 

0-668 

0-742 

0-836 

0-554 

0-643 

0-719 

0-970 

114-3 

0-645 

0-732 

0-836 

0-960 

0-448 

0-607 

0-726 

0-836 

1-049 

139-7 

0-815 

0-953 

1-113 

0-650 

0-795 

0-975 

1-090 

171-4 

0-653 

0-914 

1-143 

1-580 

0-533 

0-732 

1-110 
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TABLE  6 

TaMes  of  a/Ui  at  Traverse  PositkMis 


TABLE  6 (continued) 
Tables  of  u/Ui  at  Traverse  Positions 


TABLE  6 (continued) 

Tables  of  o/Ui  at  Traverse  Positions 


TABLE  6 (continued) 
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Fig.  2.  Sketch  of  the  working  section. 


Fig.  4.  Dimensions  of  the  plate. 


(b)  Mach  3 -6  solid  surface . 


Fig.  8.  Mach  number  distribution. 
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Fio.  19.  Tcmpcratureprofiles(M =3-6;F=0'00421). 
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Fig.  25.  Momentum  thickness  development  (M 


Fig.  26. 


Fio.27.  HvsK,(M=l-8). 
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Fig.  35.  Variation  of  dynamic  pressure  at  the  wall  (M =3-6). 
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Fig. 


Fig.38.  f7(«*)vsOog(y£4/*'».)(M=  1-8;F=000174). 
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Mixing  length  and  eddy  viscosity  distributions  (Af  =1-8;F=0-00174). 


eddy  viscosity  distributions  (M= 3-6;F=  0-00164). 
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